Z3 Quantum Criticality in a spin- 1/2 chain model 
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The stability of the magnetization m — 1/3 plateau phase of the XXZ spin- 1/2 Heisenberg chain 
with competing interactions is investigated upon switching on a staggered transverse magnetic 
field. Within a bosonization approach, it is shown that the low-energy properties of the model 
are described by an effective two-dimensional XY model in a three-fold symmetry-breaking field. 
A phase transition in the three-state Potts universality class is expected separating the m — 1/3 
plateau phase to a phase where the spins are polarized along the staggered magnetic field. The Z3 
critical properties of the transition are determined within the bosonization approach. 

PACS numbers: 75.10.Jm 



Quantum phase transitions i.e. zero-temperature tran- 
sitions driven entirely by quantum fluctuations have at- 
tracted much interest over the past decade . Quantum 
criticality can emerge from the interplay of competing or- 
ders |2J, |3| , and the existence of a quantum critical point 
offers firm ground to develop a controlled theory of the in- 
termediate regimes characterized by multiple competing 
orders Q. In one-dimensional quantum systems, a full 
description of quantum critical points can be obtained 
based on powerful analytical and numerical techniques 
such as conformal field theory (CFT) and integrability 
for the former and density matrix renormalization group 
(DMRG) and exact diagonalizations for the latter. In 
particular, the interplay between different kind of orders 
can be described, in an analytic continuum approach, 
by a CFT perturbed by several relevant terms (in the 
renormalization group sense). When acting separately, 
each perturbation yields a strong-coupling massive or- 
dered phase but the competition between them can give 
rise to a quantum critical point at intermediate coupling. 
Recently, several examples of this behavior have been 
analyzed non-perturbatively, such as the emergence of 
an Ising (or Z2) critical point in the two- frequency sine- 
Gordon model 0, 0, or in the (3 2 = An self-dual sine- 
Gordon (SDSG) model HQ. This z 2 quantum critical- 
ity has been found to be relevant for various spin chains 
Eil EH HH E3JQ_° r one-dimensional interacting 
fcrmion models jfj, lid llfi| . Richer patterns of quantum 
criticalities are expected to be realized in more complex 
situations. In this respect, U(l) Gaussian criticality and 
SU(2)fe (k = 1,2) Wess-Zumino-Novikov-Witten univer- 
sality classes have been found in the two-leg spin ladder 
with external perturbations such as a stag gered magnetic 
field 0, an explicit dimerization 

0E1LL1II3, and bi- 
quadratic interactions [21I l22j| . However, all of these crit- 
ical points can still be described in terms of free bosons 



or free Majorana fermions. An interesting question is 
whether a lattice model with only short-range interac- 
tions can yield more complex criticality. 

In this letter, we shall discuss a XXZ spin- 1/2 Heisen- 
berg chain with competing interactions in magnetic fields 
which exhibits an emerging quantum criticality in the 
three-state Potts (or Z3) universality class. We note that 
a similar critical behavior has been obtained previously in 
a numerical investigation of a quantum spin chain model 
in transverse and longitudinal magnetic fields |23j | . How- 
ever, this latter model contained a three-spin interaction 
whereas the model, considered here, involves only pair- 
wise interactions. The Hamiltonian of our model reads 
as follows: 

H = J ^2 S n ■ S n+ i + Ai ^2 Sn^n+l + ^2 ^2 SnSn+2 

n n n 



where S%, a 



x,y,z is a spin- 1/2 operator at site n 
and Ai s,J > 0. In absence of the staggered magnetic 
field h s = 0, our model corresponds to a frustrated XXZ 
Heisenberg chain with easy axis anisotropy in a magnetic 
field, the phase diagram of which has been recently ob- 
tained by the DMRG appr oach |24j . In particular, a mag- 
netization plateau [25ll2a . l27| in the magnetization curve 
of the model was obtained at 1 /3 of the full moment for 
strong enough anisotropy. This m = 1/3 plateau state 
has a three-fold degenerate ground state corresponding 
to a spontaneous breaking of the translation symmetry. 
It can be straightforwardly described in the regime of 
strong Ising anisotropy Ai^ 3> J as the formation of an 
up- up-down (TTI) structure |24|. In this letter, we shall 
investigate the fate of this to = 1/3 plateau phase upon 
application of a transverse staggered magnetic field which 
tends to order the spins antiferromagnetically along the 
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x-direction. Within a bosonization approach, we shall 
show that the competition between the up-up-down or- 
der along the z-direction and the Neel order along the 
x-direction induces a T = quantum critical point which 
falls into the Z3 universality class. The critical behavior 
of the x and z components of the spin at the transition 
will also be determined. 

The low-energy theory corresponding to the lattice 
model Q can be derived using an Abelian bosonization 
approach when Ai t 2,h g <C J. To this end, we need the 
continuum description of the spin-1/2 operators of the 
Heisenberg chain in a magnetic field in terms of a bosonic 
field $ and its dual field 9 El El El : 
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(2) 



where x — na (a being the lattice spacing) and m/2 is 
the average lattice magnetization (S%)- The Fermi mo- 
mentum hp is related to the mag netization of the chain 
as k F = tt(1 - m)/2a mHHlf. In E Q- ©> the coef- 
ficients of the expansion A\ , ^4^ are non- universal con- 
stants which can be determined numerically for a given 
lattice model j^. Within this low-energy approach, the 
m = 1/3 magnetization plateau phase of the model 
for h s = is easily recovered. In that case, the Fermi 
momentum is given by: kp = 7r/3a. The bosonizcd rep- 
resentation (J2J, yields an umklapp process in the Hamil- 
tonian JQ) stemming from the S Z S Z interacting terms, 
leading to the following bosonized Hamiltonian: 



dx 



— [ dx cosSv 7 ^^, 
■na J 



(3) 



where the Luttinger parameters v, K depend on the mi- 
croscopic parameters of the lattice Hamiltonian with 
K decreasing as A 1.2 increase. The low-energy effective 
Hamiltonian © thus corresponds to the f3 2 — 18nK sine- 
Gordon model. The sine-Gordon form of Eq. J3J can also 
be justified on symmetry grounds. Indeed, translation by 
one lattice site T a is described by: $ — > $ + V2~rr/3 in 
the bosonization approach j^, and cosS-v/^-*! 5 is, among 
the operators left invariant by this symmetry, the one 
with the lowest scaling dimension. The field theory de- 
fined by Eq. |J3J becomes massive when K < 4/9. In 
that regime, which can be reached provided Ai, A 2 > 
are sufficiently large |29j, the bosonic field $ is pinned 
on one of the minima of the sine-Gordon model © i.e. 
($) = pV2tt/3 for A > 0, p being integer, so that a mag- 
netization plateau at m = 1/3 of the saturation magne- 
tization is formed. Taking into account the compactified 



nature of the bosonic field <1> (<& ~ <f> + \/2w) in Eq. J5J , 
we observe that there are three non-equivalent ground 
states: ($) = 0, y/2n/3, 2^/2~n/3 that are transformed 
into each other by a translation of one lattice site. This 
three-fold ground-state degeneracy is the thus the result 
of the spontaneous tripling of the unit cell in the plateau 
phase, as can also be seen from the decomposition |J5J, 
which shows that (S*) oscillates with a period 3 while 
(S%) = in the plateau phase. 

The low-lying excitations of this m = 1/3 plateau 
phase are massive kinks and antikinks interpolating be- 
tween consecutive ground states, i.e. configurations of 
the field such that $(+00) - $(-00) = y/2n/3. 

These massive kink excitations thus carry fractional S z 
quantum number: S z = J^^dx d x ^/V2w = 1/3. 
All these properties are in agreement with the Ising 
anisotropic limit: Ai,A2 3> J,h s . In that limit, the 
ground state is three-fold degenerate for a magnetiza- 
tion per site (S z ) = 1/6 with an up up down (TTI) 
structure invariant by a translation of three lattice sites. 
Moreover, the low- lying excitations of this structure are 
known to be massive domain walls carrying a fractional 
S z = 1/3 quantum number |24l l30j in full agreement 
with the bosonization approach. Up to now, we have 
been considering the case of a field tuned especially to 
the value H = H c so that the state with m = 1/3 is the 
ground state. If the field starts to deviate from the spe- 
cial value H = H c , it will generate a chemical potential 
for the kinks and the antikinks (H — H c )d x <&. However, 
if this chemical potential is smaller than the gap to cre- 
ate kink excitations, no kink will appear in the ground 
state, and the magnetization will remain at m = 1/3. For 
larger deviation of the magnetic field, a commensurate- 
incommensurate transition will occur and the magnetiza- 
tion will deviate from m = 1/3. This is the origin of the 
presence of the magnetization plateau for a finite interval 
of the magnetic field. 

Let us now investigate the stability of the m = 1/3 
plateau phase upon switching on the transverse staggered 
magnetic field h s . Clearly, for h s —> 00 the spins will be 
polarized along the x-direction and the m — 1/3 mag- 
netization plateau will be destroyed. The open question 
concerns the nature of the phase transition between these 
two strong coupling massive ordered phases. We begin 
with the case H = H c , and later on turn to H ^ H c . Let 
us write down the bosonized Hamiltonian corresponding 
to the model JD when h s <C J: 

n = V -J dx [(d x e) 2 + (d x $f] (4) 

• J dx cos ^2-k/KQ - — - J dx cos VIEttK^, 

where we have rescaled the bosonic fields by the Lut- 
tinger parameter K. The model Q is a generalization 
of the sine-Gordon model © obtained by adding a ver- 
tex operator which depends on the dual field 0. The 
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scaling dimensions of the two perturbations in Eq. @ 
are respectively A hs = 1/(2K) and A A = 9K/2. For 
1/4 < K < 4/9, a regime which can be reached owing to 
the presence of two independent couplings A 12 > in 
Eq. QJ, the two interacting terms are strongly relevant 
perturbations so that a phase transition is expected. 

Let us first give an argument for the existence of a 
phase transition for finite and nonzero h s . For h s = 0, 
i.e. the m = 1/3 plateau phase, the model (@J reduces to 
the pi 2 — 18ivK sine- Gordon model on the field 4> with 
massive solitons as elementary excitations. The soliton- 
creating operators in this sine-Gordon model with topo- 
logical charge m read as follows: 



O m = 



cxp I I — 



(5) 



which, in the path integral representation, creates a dis- 
continuity of the field <I> ($ — > $ + my / 2n/K/3) along 
a contour of integration ( "Dirac string" ) which enters in 
the definition of the dual field [3l| . Upon switching on 
a non-zero value of h s , we notice that the perturbation 
cos{y/2-K / K&) in Eq. Q is locally mutual with the soli- 
tons © of the f3 2 — 187T\ftf sine-Gordon model on the 
bosonic field $ and can be interpreted as creating or an- 
nihilating a packet of three elementary solitons. The lo- 
cality of the interaction term cos(y / 27r / KQ) with respect 
to the solitons JBJ of the sine-Gordon model © implies 
that no phase transition occurs as an infinitesimal stag- 
gered magnetic field is switched on. The m = 1/3 plateau 
phase is thus stable and extends to a non- vanishing value 
of h s . Let us now make similar considerations for the op- 
posite regime (h s ^> A) of the model J3J where we are now 
perturbing around the 1 = 2tt/K sine-Gordon model on 
the dual field ©. The ground state of the model (@J with 
A = is characterized by (5*) = and ((-l) n S%) ^ 
as it can be easily seen from the decomposition @. In 
the same way as for the m = 1/3 plateau phase, this 
Neel ordered phase along the x-direction is stable up to 
a strictly positive value of A. The conclusions for the 
two limits considered here are only compatible if a quan- 
tum phase transition occurs at a finite non-zero value of 
the staggered magnetic field h* between these two strong 
coupling massive phases. 

The nature of the emerging quantum criticality at h* 
can be elucidated by noting that the low-energy field 
theory (0} also describes the critical properties of the 
two-dimensional classical XY model with a three-fold 
symmetry-breaking field defined by the following lattice 
Hamiltonian 1321: 



Hz 3 =-J ^2 cos(ip r - <p T i) - /i^cos(3 <p r ) , (6) 



(p T being the angle of the unit-length rotor at site r of 
a square lattice, and h breaks the continuous 0(2) sym- 
metry of the XY model down to a discrete one, Z3. For 



h — > +00, the model Jfjjl reduces to a Z 3 clock model 
which is equivalent to the three-state Potts model. More 
generally, since at zero temperature the term h cos 3</j se- 
lects ip = — 2tt/3, 0, 27r/3 as the possible ground states, 
the model © is expected to display a phase transition as 
a function of temperature in the three-state Potts (Z3) 
universality class. This Z3 criticality has indeed been 
confirmed numerically by means of series analysis |33| . 
and exact diagonalizations on finite samples [34|. It has 
been known for some time |35j that the universal prop- 
erties of the model (@ are captured by the Hamiltonian 
0. In particular, the first two terms in Eq. (0J consti- 
tute the continuum description of the 0(2) XY model, 
with the cosine of the dual field in Eq. J3J creating and 
annihilating topological vortices of the field cp r . In addi- 
tion, the vertex operator with coupling constant A in Eq. 
(0J represents the continuum limit of the three-fold sym- 
metry breaking perturbation. The phase transition of the 
model at /»*, separating the m = 1/3 plateau phase 
from the Neel ordering along the x-direction, should thus 
belongs to the Z3 universality class. 

Some insights of the physical properties of the Z3 qua- 
tum critical point might be gained by considering an ana- 
lytic approach a la Luther-Emery |3jj of the generalized 
sine- Gordon model (J2J. At the special value K = 1/3 
and by fine-tuning the coupling constants, the model (@J 
simplifies as follows 



Hsdsg 



dx 



(d x $>) 2 + (d x &) 2 



J dx cos (V67T + cos (V6w ,(7) 



which exhibits a self-dual $ <-> O symmetry and will be 
referred in the following as the (3 2 = 6ir SDSG model. 
It has been shown recently in a non-perturbative fashion 
that the (3 2 = 6tt SDSG model displays an emerging 
Z3 quantum criticality 0, • The universal properties 
of the model (0 can be determined using the ultraviolet- 
infrared (UV-IR) transmutation of some fields of the 
massless flow of the 1 = 6-7T SDSG model which has been 
obtained from a Toulouse limit analysis [|J. In particu- 
lar, it has been shown in Ref. IJthat the vertex operators 
exp(±i v /27r/3 $) of the (3 2 = 6ir SDSG model identify 
in the IR limit to the two Z3 spin fields a\ and 02 = <r\ 
with scali ng d imension 2/15 and q = ±1 Z3 charge re- 
spectively |38| . Using the continuum representation J3J), 
we then predict the leading asymptotics at the Z3 critical 
point (h s = h*) of the spin-spin correlation function in 
the z-direction: 



({S* (x)S z (y))) 



([S z (x) - m/2a] [S z (y) - m/2a}) 
cos (2n (x — y) /3a) 



y\ 4/ 



15 



(8) 



In the same way, it is possible to show that the opera- 
tor cos(\/67rO) transmutes in the IR limit to the energy 
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operator e of the three-state Potts model with scaling 
dimension 4/5 [3^. We thus obtain from Eq. (J2J the 
leading asymptotics of the spin-spin correlation function 
along the x-direction at the Z 3 critical point as: 



«S* (x)S*(y))) 



\x~y\ s / 5 ' 



(9) 



The spin-spin correlation functions I|8I9|) thus display an 
algebraic behavior with universal exponents characteris- 
tic of the Z 3 universality class. When the system is tuned 
away from the self-dual point, we have to consider which 
relevant operators can be generated. Both the operators 
(d x &) 2 and cosa/67tG — cosv / 67r ( I > are mapped onto the 
Z3 thermal operator e in the IR limit. Can another rel- 
evant operator be generated besides the thermal one ? 
The entire content on the Z3 Potts model is known [ssj 
and the thermal operator is the only relevant operator 
having zero conformal spin and preserving the global Z3 
symmetry of the model. The other relevant operators of 
the Z3 Potts model carry a non-zero conformal spin and 
only play a role when one deviates from the center of the 
m = 1/3 plateau phase as we will discuss below. As a 
result, all deviations from the self-dual point reduce here 
to tuning the temperature in the Potts model away from 
the critical temperature. The duality mapping <i> «-> 0, 
K *-* 1/K can therefore be seen as the duality transfor- 
mation of the Potts model between the high and the low 
temperature phase. The presence of the two indepen- 
dent coupling constants A2 and h s of Eq. |Q should be 
enough to reach the Z3 critical point. 

So far, we have only discussed the case of a mag- 
netic field H = H c . If we now modify the mag- 
netic field to deviate slightly from the center of the 
m =1/3 plateau phase, a Zeeman term Hz — —{H — 
H c )(K/27r) 1 / 2 Jdx d x $ should be added to the low- 
energy Hamiltonian The resulting Hamiltonian de- 
scribes a two-dimensional asymmetric XY model in a 
three-fold breaking field: 

Hs = -jy^cos(y r -y r+e „-^)~ fe y^ cos ( S( Pr) , (10) 

r,/i r 

with 8^ ~ (H — H c )e x ■ e M , (/1 = x, y). When h — > 00, 
this model reduces to the asymmetric three-state clock 
model which displays commensurate, paramagnetic, 
and incommensurate phases |39T.l4oj| . Close to the Z3 crit- 
ical point, the low-energy effective Hamiltonian density 
corresponding to Eq. fTT))) is given by |4lj : 

H 5 = Hz 3 + £$(2/6,7/5) + 5*^(7/5,2/5) + (h s ~ K) 6, (11) 

where $(2/5,7/5) and $(7/5,2/5) are Z3 operators with scal- 
ing dimension 9/5 and conformal spin ±1. The model 
f| 1 1 1> with 5* = and h s — h* is known to display critical 
behavior in the chiral three-state Potts universality class 
due to the presence of a non-Lorentz invariant pertur- 
bation 4.1.1. Such a chiral transition has been obtained 



recently in the phase diagram of a one-dimensional bo- 
son model with competing density- wave orders |42| . For 
<5 = S* , there is a possibility of having a Lifshitz point be- 
tween a chiral Z3 and a standard Z3 critical lines. In this 
respect, it will be interesting to investigate numerically 
the phase diagram of the lattice model (JIJ to further shed 
light on this issue. 
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